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ABSTRACT 


High fidelity elastic system modeling algorithms are discussed. The particular system 
studied is the Space Shuttle Remote Manipulator System (RMS) undergoing full articu- 
lated motion. The model incorporates flexibility via a methodology the author has been 
developing. The technique is based in variational principles, so rigorous boundary condi- 
tion generation and weak formulations for the associated partial differential equations are 
realized, yet the analyst need not integrate by parts. The methodology is formulated usine 
vector-dyad notation with minimal use of tensor notation, therefore the technique is believed 
to be affable to practicing engineers. The objectives of this work are to: 

1. Determine the efficacy of the modeling method. 

2. Determine if the method affords an analyst advantages in the overall modeling and 
simulation task. 

Generated out of necessity were Mathcmatica algorithms that quasi-automate the modeling 
procedure and simulation development. The project was divided into sections as follows: 

1. Model development of a simplified manipulator. 

2. Model development for the full-freedom RMS including a flexible movable base on a six 
degree of freedom orbiter. A rigid-body is attached to the manipulator end-effector. 

3. Simulation development for item 2. 

4. Comparison to the currently used model of the flexible RMS in the Structures and 
Mechanics Division of NASA JSC. 

At the time of the writing of this report, items 3 and 4 above were not complete. 



INTRODUCTION 


Material bodies are inherently of a distributed mass and elasticity nature. Analyst have 
realized this fact since the early days and developed tools to model these distributed ef- 
fects [1], Engineers, challenged with the task of making devices work in a reliable, energy 
efficient, and inexpensive manner, have been gradually increasing the fidelity of their models 
by incorporating the distributed properties. The ability to study these high fidelity models 
grows with the increasing computational capabilities of inexpensive computers. 

The literature is teeming with ever-improving ways to model the distributed effects [1], 
There are a diverse cross-section of techniques. Some are intuitive to a design engineer [2. 
3, 4, 5], while others are mathematically elegant but beyond the training of many practicing 
engineers [6, 7]. The purpose of this study is to examine the efficacy of the author's attempt 
at developing a rigorous yet usable method for modeling complicated systems [5] . 

METHODOLOGY 


Present Capabilities 

Based on discussions. 1 the author understands that the fidelity of the model for the 
present Shuttle Remote Manipulator System (RMS) simulation is limited to small amplitude 
vibrations about any “snap shot” configuration of the system. This limitation manifests 
itself because of the linear finite element (NASTRAN) model used as the progenitor for the 
modal basis functions. Therefore, RMS slewing maneuver studies are not within the fidelity 
of the linear model. There exist techniques which allow an analyst to study the slewing 
maneuvers of systems like the RMS, but these modeling techniques are computationally 
expensive and/or hard to understand [1], therefore they are not always implemented by 
practicing engineers. The author believes the technique discussed below gives analysts a 
familiar yet powerful modeling tool. 

New Capabilities 

The main motivating factor for the development of another modeling method was the need 
to easily derive complete models of complex elastic systems [1. 4, 8, 9, 10, 11, 12]. Although 
the method discussed herein is still relatively mathematically intense (compared to an equal 
number of rigid bodies), it has a predisposition for symbolic manipulation. Another impetus 
for this work is that a simple method may make it possible to bring rigorous flexible svstem 
modeling out of the academic domain and into use by product designers. Another catalvst 
for this effort is that a simple (ultimately an automated) method will make it possible for 
researchers to rapidly regenerate models based on new continuum assumptions. 

The approach is variational in nature. It retains most of the attributes of the analytical 
approach (i.e. Hamilton’s principle), but eliminates most of the pitfalls, such as the need to 
use Lagrange multipliers for constraints, and excessive algebra. The methodology is vector 


‘Orientation meetings with various engineers from the Structures and Mechanics Division of JSC. 




Figure 1: Two Link Flexible Manipulator 


based and requires the analyst to perform operations comparable to the operations required 
for implementing Lagrange’s equations. However, it is claimed that the net algebra with 
the method herein will be less than the net algebra associated to Hamilton’s principle or 
Lagrange s equations. Analysts familiar with Kane s [13] form of d’Alembert’s principle will 
find the technique affable. The complete derivation of the method is shown elsewhere [1, 5], 

SIMPLIFIED MANIPULATOR MODEL 

Preamble 

In this section the equations of motion for a flexible two-link planar manipulator will 
be derived. This section is included for demonstrative purposes. The procedure that is 
used on the full RMS model is used on a somewhat simpler model so the reader can follow 
the steps involved. This simpler system was chosen because: a) it is non-trivial due to its 
distributed elasticity, b) its planar nature allows for heuristic equation verification, and c) 
Hamilton’s principle can be readily applied to it. This example will demonstrate some of the 
qualities of the new methodology, such as: a) its systematic nature, b) its resulting closed 
form equations, and c) automatic boundary condition generation. 

The system is shown in figure 1. The domain of each beam is one dimensional. The 
independent coordinates are x n and measured from the root of beam £, and B,. re- 
spectively, along the undeformed neutral axis of each beam. The “special” point of beam 
& i labeled a, (the pivot) and s 2 (the joint between the two beams) is for beam B 2 . The 
coordinate frames, denoted with 5, and B 2 , are attached as shown in figure 1. At the root 
of each beam (£?i and B 2 ) there are massless hubs to which torques .V/, and \l 2 are applied. 
The angular position of frames B\ and B 2 are <jq and q$. Beam deflections are measured 
with u u (£ u ,t)b n and u 2 i(x 2U t)b- n (elongation), and with u X2 {x u ,t)'b r2 and u 22 ( x 2 i . t)b 22 
(flexure), as shown in figure 1. The beams have mass per unit length p, total lengths are 
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1 1 aud h C f r ? SS se ^ ti ° ,lal . area ' 4 . area moment of inertia I, and Young’s modulus E It is 
assumed that large deflections and rotational inertia is pertinent, but not shear deformation 

the be r s are mod i ed with Rayieigh ^ c^^tiousTeS; 

beam are assumed symmetric about the neutral axis. The intrinsic mass moment of inertia 
of the cross-section will be taken as I lo = pll ,f 3 . * 3 is normal to the plane of the problem 


Kinematics 

The geometry of the motion for the system is as follows. The 


are: 


position vectors of interest 


°f*« = 0 

1 1 = ' r = (-Tn + iiii ) 6 j] + u, 

° rJJ = (^i + q' 2 ) b u + q' 3 b u 

05 = ( x 21 + 7 * 21 ) ^21 + U 22 6 - 2 .) 

r ~ (^1 "b (]■}) 611 + <73^12 + (x 2 i + U'2l ) 6-21 

+ “ 22^22 


( 1 ) 

( 2 ) 

(3) 

(4) 

(0) 


The angular velocity of frame B, and B, and intermediate frames (in the cross-section) /, 
and I2 are: 71 


jV-> b, ; 

^ = Ui 6| 3 

Kr ,* a - 


u7 2 — ( U 1 + uj, + 125) 613 


JV- 1 , _ / . ^ 2 Ui 2 \ ; 

* - r + 5ns u 


13 


^ /j = ( «1 + U' + W 5 + 


dx^dt 


bn 


( 6 ) 

(7) 

(«) 

(9) 


The generali2ed and pseudo-generalized (denoted with a /) coordinates and veloci 
denned as: 


U] = 91 , u 5 = q 5 

<12 = 5ll(£l,0, <7^ = 5,2(i 1 ,0, <?' = - 12(Ll ^ } 

ax n 


ities are 
( 10 ) 

(u) 


71 1 1 ( Z. 1 , 7 ) 

dt 

= A t _ 

d.r u dt 


u 2 = <h = 


u 4 = q 4 = 


, U3 = </., 


_ dtIi 2 (Z.|, t) 
dt. 


( 12 ) 
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The absolute velocity of the point .Si and s 2 are also required. Since the system is rotating 
about .s i : 

°vj'=0 (l.T) 

and 

°Vtf 2 = (u£ u 1 7.3 ) &n + ( u 3 + u \ (^i + 92 )) Wi (14) 

where the pseudo-coordinates and pseudo-speeds have been defined above. 

The absolute acceleration of the differential beam elements for beam B\ and B 2 can be 
written as: 


and 


O- lo , _ f d 2 Un . ~ duu , , ~ , A ; 

a A — I ( X ll 4" ^ll) U\ 1 V\\ 

( d 2 u u _ du n „ A- 

+ 1 -QP f (Xu + Uji j Ui + 2 U] — ^ Ui/Uj I bn 


15) 


5 a^ = («2 - - 2 u,«3 -(/,,+ <7-2) u?) £>n 

+ (U 3 -I- (Ii + < 7 - 2 ) ui + 2 ujU 2 - (fouf) b u 

+ ^ ~ (ui + U 4 + u 5 ) u 2 - 2 - 


du- 


2 (ui + u' 4 + us) — — (a.' 2 i + u-j 1) (uj + u 4 + u 5 )' ) 621 


8t 

' d 2 u n 

+ I ' 7T7T - + (X 21 + U- 2 l) (^1 + W 4 + U 5 ) 


dt 2 


du 21 


+2(Ui + u\ + Us) — ^ U 22 («l + 


u\ + u. 5 ) 2 ^ 


br. 


■22 


16) 


respectively. For this problem z oe for each beam is at the centroid of the cross-section. 

The methodology also requires the calculations of the "preferred directions” for the vari- 
ations (pseudo and ordinary), namely the partial speeds. They are determined by inspection 
of the velocity equations and given in table 1. The partial velocities for the field equations 
are given as: 


d Bl u ,x 

— lu > 

d B2 u h 

du^iM 


dun,\ t 

d 3 x v B \ x 

= b u . 

d sx v B i' 

d u \\ ,t 

du | 2 c 

d**v B ? 

= L. 

d 32 v p ? 
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°Vjt' 

°r 


- V J 62 

3 
du \ 

0 

~ < 73^11 + (^i + q' 2 ) ~bi 2 + (iiibn 

bi 3 

b'2:i 

3 

3uif 

0 

"T 11 

0 

0 

3 

0 

h-i 

0 

0 

a 

3u\ 

0 

0 

0 

^23 

3 

3 

0 

0 

0 

b'2'.i 


Table 1 : Partial Velocities for Pseudo and Regular Coordinates 


The strain energy density functions for the beams B x (i = 1,2) are (assuming large 
deflections): 




n d °^‘ f v i 

0 = ‘ i Ffi ' " /b ‘ 


du i 
du i 


+ 




[Ts, — Jb 1 


du i 

Fb 2 — i * B 2 ] H ^ — * [ Tb 2 ~ Jb 2 


for u l? and 


0 = ^ • [fg, - / Bl ] + . [fg, - jg,] 


du-, 
d°vf 
du 5 


Cs 2 - I b 2 H 


du 5 

d' v jj B 2 




[7 b 2 - Jfl., 


(IS) 


The torques applied to the massless hubs are: 

T\ — M\ biz and Ti = M^by 3 ( 1 9 ) 

on Bi and B 2 , respectively. The equations of motion can now be written down. 

Ordinary Differential Equations 

The ordinary differential equations, governing the angular positions are: 


( 20 ) 


( 21 ) 


for u 5 . The partial velocities are defined in table 1 . The forces and torques (applied and 
inertia) are defined as: 


A = j ‘ { KF lie + v r F df ) dn t 
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+ / fa' The + X f F de ) da f 

Jan. v ' 

= f p°dx u dQ e 

Jn. 

= J ['*r‘~ x faj he + vj de ) 


+'H e Tk e 4- X e T de rffi 


/ K r ’~ x fa e F he + V e F dl 
Jan. L ' 


+U.T ke + T> e f de 1 (fcr, 


■h = / ,e f' ‘ oe X ft'ay'* + 

•/n* 


‘^r x p f °aj°' + • A a '* + ' V J '• x /,„ • 


^ /* . a*- /, 


The final form of the differential equation are found by taking the indicated dot products 
and are not displayed here. 


Partial Differential Equations 

The field equations governing elongation and bending for B u 

_ d ( dVi \ . . 


dx n 2 \dx 


_P UlUli _ “ Ul_ at ( Xu + Ull > 


for elongation, and: 


d ( dv x \ d 2 f dv. 


dx u V^^i2.t/ ^11 \ dtij2 




ea [ ^in + i (*hi 


J>L ( 

dx 2 u \ dxj x 


2-8 



. f ^*^12 , 0 du l2 _ __ 

\ ^ ( X M + U II ) “l - u u uf 

d 17. d 3 5i2 \ 

“1 + ~T I pj 


dx 


11 


+pl u, + 


dx u dt 2 y 

&% 2 x 


<?x„dJ (^3 • *13)(&13 X in) . 6,3 


for bending. 

O 

At x u = 0, the boundary conditions for B x 


are: 


u \ i = u V2 = u , 2 | = 0 

The boundary conditions at z„ = L x for B x are: 


^ , ( du u 1 / dui x 2> 


^ u i i.i 


= EA 


dx u + 2 {dx n l ) 5,1 


‘12 


for elongation, with: 


m 


(24) 


(25) 


for shear, and 


9u = *13 • [l\ ■ '7* '■ + *C3 '« x ■ "u '<) 
, dV\ d ( d\\ \ 


^“i2,i dx n \du 


‘12,11 , 


_ . d 3 u u \ 

~ I U 1 + T-T I pi 


dx u dt 2 


d 2 ui. N 2 


+pI + d^k ) (bv) ' ‘ l3 )&3 X 2 »3) • *13 


+ EA ( + A ( *\ du 


12 


dx 


11 ^ \ 3 £ 11 


dx 


11 


d 


dx 


El 


ii 


d^U \ 2 
0*11 , 


P- = ^ 

^ U 12.11 dxi, 12 


11 


defined^* m ° ment ' Tll? illtrillsic terms of the boundary conditions at r 


11 = 


(26) 

(27) 
L x are 


' _ ggg fa f 1 r- 

du\ • r*' - + -^r- • [ r » 


d x u) B ' 


Jb> 
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d°v& 

dll'; 


r — — 1 

d y Co B > 

- 1 

Fbi - Ib 2 

+ c)u' 2 

Tg 2 — Jb 2 


and 


9 12 = 


d%? 

du-x 


r - - n d J ^'~ e ' 

[^. “ ^ ' (^Bi “ '•'Hi 


I fr T 1 , **&* [f T 

+ ~d^~ ■ e B ’ ' ,B >] + ' l 7 ®’ “ Jb ■■ 


and 








du 4 


The field equations for the second member ( B 2 ) are: 


0 = 


d ( d V 2 


dx 2 i \$t * 21 j 

d 


- p°a ■ 621 


c)x 


21 


£/i 1 |Hhi + i 


()x 


21 


2 V dx 


21 . 


-p (u 2 — u\(j :} — 2 uja 3 — u, (L\ + #/-2 ) j (&11 • &21) 

-p [1X3 + tii (£1 + 72 ) + 2 uiu' 2 — Uj 7 jj (612 ■ £>21) 
^ 2 U 2 l 


-p 


dt 2 


— (tii + U4 + tis) U22 


-2 (ui + u\ + u 5 ) - (ui + u\ + u 5 ) 2 (x 2 i + U21) 


for elongation, and: 


0 = 


d Vi 


d 2 


d Vi 


dx 2l \du. 2U J dx 2 n \du Tt 


22,11 


P°0 !> • b> 2 


d 


dx 


21 


0 


Dx 


hi ' (X • A a /2 + jV ^ /2 x 


^4 [ + 1 ( d “*>Y\ du-n 


21 


cfj- 


21 


2 V C/X ; 


21 


Ox 


21 


( 28 ) 


( 29 ) 


( 30 ) 


( 31 ) 
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d 2 


El 


dx%i V dz 2 > 


21 


p h U] q 3 2 uiu' 3 - U 2 (h + q' 2 )j (& n • b 2i ) 

~P [ u 3 + u i (^1 + q 2 ) + 2 u\u 2 - u\q 3 j (b u • 6 22 ) 
d 2 u i2 / , ; 

^2 (“• + ^4 + “5) (*f 21 + u 2 i) 


<40, 


+2 (u, + u\ + u 5 ) - (u, + u\ -f u 5 ) 2 u 2 . 


+ - 


dx 


21 


' ■ , •/ . ■ C? 3 U22 \ 

, Ul+ “' + “ s + ^? k 1 


+ pz lu \ + u 4 + u$ + 7 


3 x 2 \ 3t 


[ bz : 3 ' 4 2:3) (^23 X *23) ’ ^23 


for bending. 

The boundary conditions for 5 2 at x 2 t = 0 are: 

U21 = U22 = 222,1 = 0 

The boundary conditions at x 2] = l 2 are: 


0 = EA 


du^ 1 / du 2 2 
dx>, + 2 \ ~dx 


21 . 


for elongation, and: 


0 -(*+*; + *. + £*)* 


+ P Z + u' + u 5 + j-!L \ { l i3 . l 23 ){ l 23 x - 23) . ii3 


+ E 4 | - Un + I ( _ 2 \ 11 

dx n 2 ^x 21 / / tfx 2I dx, 1 ' ' 


21 ox 2 , 


dx 2 n 


for shear, and: 


0= El 


dx 2 n 


(T2) 


(T3) 

(34) 

(33) 


(3(j) 


(37) 


for moment. 

Initial conditions for q u q s , u 5 . 5 n , u 12 , u ai . u 2 2, %*-, and % must aIso 

be speched. The kinematic differential equations for q , and * are given in equation 10. 
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Discussion 


As can be seen by the presentation above, the methodology facilitates the process of 
writing equations of motion for complex systems. One can see that the rigorous natural 
boundary conditions that are generated via variational principles are present but the analyst 
does not have to integrate by parts. For comparison of the technique with Hamilton's 
principle see [14] relative to the problem above. For more complex systems, such as those 
exhibiting nonholotiomic constraints, see [1, 15]. For systems undergoing contact/impact, in 
a hybrid parameter fashion, see [I, 16. 17], 

Apparently some of the lingering questions are whether or not the method is attractive 
to practicing engineers and can the symbolic form of the equations be put in a form suitable 
for simulation. The author claims yes to the later question and only time will tell on the 
former question. Massaging the equations into a simulation are the topic of the next sectiou. 

With regards to getting numbers from the equations, it is possible to put the equations 
in a weak form so that the complicated boundary conditions are absorbed into an integrated 
form of the partial differential equations. Then one has only to choose an appropriate 
function for the test function used to cast the problem in its weak form. For problems like the 
manipulator above, Rayleigh-Ritz discretization is probably sufficient. For more complicated 
continuum bodies, finite element discretization is probably appropriate provided the analvst 
work with the problem in its weak form so the appropriate boundary conditions are included. 

RMS MODEL 


Model Description 

The main subject of this report is the application of the method described above to 
the Space Shuttle RMS. The RMS is modeled as a system of rigid and continuously elastic 
bodies, a hybrid parameter mechanical system. The system is broken down as follows (see 
figure 2). 

The orbiter is taken as a six degree of freedom rigid body. The RMS base (assumed to 
be rigid and labeled B ) is attached to the orbiter via small displacement small angle springs 
which approximate the elastic nature of the orbiter. The RMS shoulder yaw motor and 
housing are assumed to be a rigid body (body Sy). Connected to Sy via a nonlinear spring 
and motor control action is the RMS shoulder pitch body Sp with its actuator. Attached 
to Sp is the first elastic boom. This boom, as with all booms in this model, is assumed to 
be a continuum in which y and r deflections along with axial rotation are modeled. The 
beam model is a Rayleigh beam (intrinsic cross-sectional inertia incorporated) with small 
deflections. The next body in the chain is the elbow pitch motor and housing. It is attached 
through a nonlinear actuator to rigid body EL. Attached to EL is the next flexible boom, 
modeled like the first boom. Attached to the second boom is the wrist pitch motor and 
housing. It is attached to body P. Boom three is anchored by body P. Boom three is 
elastic and modeled as described before. At the tip of boom 3 is the motor and housing for 
the wrist yaw action. The wrist yaw body is rigid and labeled ) Attached to the yaw budy 
through a nonlinear actuator is the wrist roll motor and housing. This is followed by the 
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fourth elastic boom (the end effector). Attached to the end effector is a rigid payload. 

Closed Form Model 

In order to attempt writing a closed form model for a system as complicated as the RMS 
described above, a symbolic manipulating assistant is desired. The author has access to 
Mathematica so this is the tool that was utilized. Before getting into the description of the 
algorithms developed, justification for the effort will be outlined. 

Why should an analyst develop closed form models when there exist other tools that 
seem to adequately model these systems? The author believes that using tools that are 
traditionally from the structural analysis realm such as NASTRAN models unnecessarily 
limit the model to the linear motion about some configuration. It is felt that if the approach 
of writing complete models first (then reducing to linear if desired) is feasible, in a timelv 
manner, then engineers will utilize these more exact models. In order to facilitate the clock, 
computer aided modeling is desired; Mathematica is an excellent tool for this process. An- 
other advantage to working directly with the closed form model is that the “zero times zero” 
multiplications that arise in straight out matrix models are avoided. Also repetitive multi- 
plications and additions are readily recognized and can be assigned to a memory location for 
instant recall. This tight code will make running these complicated models more feasible. 

Mathematica Algorithms 

Mathematica algorithms were developed to mimic the procedure outlined in the previously 
discussed simplified model of the RMS. The standard notation for Mathematica was adjusted 
so as to mimic engineering vector notation. Then algorithms were developed that recognize 
the vector dot and cross products, the triple products, and other identities. Differentiation 
of vectors in multiple coordinate frames was defined. Standard order for the symbols was 
defined so symbolic cancelation was facilated. Function that aid in the gathering of terms, 
the distribution ot terms, and general manipulation were developed. At this point these 
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algorithms are used via a Mathematica notebook running on a NeXT computer. They are 
not limited to this computer system because the notebooks are portable across multiple 
computer systems. An example of how one enters symbols for manipulation is shown in the 
appendix. 

RMS Model Status 

Presently the modeling procedure is not complete. All stages ul the development are 
complete up to the point where the actual differential equations suitable for output to FOR- 
TRAN format are formed. All the appropriate d'Alembert, forces and torques have been 
calculated along with the appropriate partial velocities, and the weak formulation. Unfortu- 
nately the approach taken thus far is very memory hungry so the workstation is using a lot 
of virtual memory which is time consuming. Refinements to the procedure and algorithms 
are made in real time and the memory and time consumption problems are being reduced. 

Model Shake Down 

Comparisons of the aforementioned model with existing models will be made upon com- 
pletion of the modeling procedure. It is intended that the efficacy of the technique and model 
will be tested via metrics such as accuracy, calculation speed, and generality. 

SUMMARY 

The rudimental aspects of a procedure to rigorously model complicated systems in a 
timely manner have been developed. The modeling technique is based on a variational 
principle based approach for writing the equations of motion, augmented with computer 
aided modeling algorithms written in Mathematica code. The tools developed are being 
applied to a complicated RMS model in order to establish the efficacy of the modeling 
technique. The technique shows promise because of its rigor, but the details of the computer 
aided algorithms need refinement. Numerical studies have yet to be performed. 
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Appendix.ma 


1 


This is the symbol manipulations for the robotic arm flex model. 


Angular Velocity and Acceleration, omega (W ) are the generalized speed. 


Newtonian Frame (N) 

Orbiter Frame (O) 

o [1] : =unitVector [O, o, 1] 
o [2] : =unitVector [O, o, 2] 
o [3] : =unit Vector [O, o, 3 ] 

NwO=omega [N,0] =Wol o[l] + Wo2 o[2] + Wo3 o[3] 


Wol o + Wo2 o + Wo3 o 
12 3 

NaO=DvDt [N, NwO] //Simplify 


Wol o +■ Wo2 o + Wo3 o 
12 3 

Manipulator Base Frame (B) 
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2 


b[l] : =unitVector [B, b, 1] 
b[2] : =unitVector [B, b, 2] 
b[3] : =unitVector [B, b, 3] 

OwB= omega [0, B] =Wbl b[l] + Wb2 b[2] + Wb3 b[3] 

AAA 

Wbl b + Wb2 b + Wb3 b 
12 3 

NwB= omega [N, B] =NwO + OwB 


Wbl b + Wb2 b + Wb3 b + Wol o + Wo2 o + Wo3 o 
12 3 12 3 

NaB=NaO + DvDt [ B , OwB ] + NwB >< OwB 


Wbl Wol o >< b + Wb2 Wol o >< b + Wb3 Wol o >< b + 
11 12 13 


Wbl Wo2 o >< b + Wb2 Wo2 o >< b + Wb3 Wo2 o >< b + 
2 1 2 2 2 3 


Wbl Wo3 o >< b -i- Wb2 Wo3 o >< b + Wb3 Wo3 o ><. b + 
3 1 3 2 3 3 


Wbl b + Wb2 b + Wb3 b + Wol o + Wo2 o + Wo3 o 
12 3 12 3 

zl=Coeff icient [NaB, o [1] ><b [1] ] 
z2=Coef f icient [NaB, o [1] ><b [2] ] 
z3=Coef f icient [NaB, o[l]xb[3]] 
z4=Coef f icient [NaB, o [2] ><b [1] ] 
z5=Coeff icient [NaB, o [2] ><b [2] ] 
z6=Coeff icient [NaB, o [2] ><b[3] ] 
z7=Coeff icient [NaB, o [3] ><b [1] ] 
z8=Coef f icient [NaB, o [3] ><b [2] ] 
z 9 =Coef f icient [NaB, o[3]xb[3] ] 


Wbl 

Wol 

Wb2 

Wol 

Wb3 

Wol 

Wbl 

Wo 2 

Wb2 

Wo 2 

Wb3 

Wo2 
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